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Abstract 

In this paper, we study the probabiUty that a dense network confined within a given geometry is fully connected. 
We employ a cluster expansion approach often used in statistical physics to analyze the effects that the boundaries of 
the geometry have on connectivity. To maximize practicality and applicability, we adopt four important point-to-point 
link models based on outage probability in our analysis: single-input single-output (SISO), single-input multiple- 
I output (SIMO), multiple-input single-output (MISO), and multiple-input multiple-output (MIMO). Furthermore, we 

derive diversity and power scaling laws that dictate how boundary effects can be mitigated (to leading order) in 
confined dense networks for each of these models. Finally, in order to demonstrate the versatility of our theory, we 
^3 analyze boundary effects for dense networks comprising MIMO point-to-point links confined within a right prism, 

^ ^ a polyhedron that accurately models many geometries that can be found in practice. We provide numerical results 

for this example, which verify our analytical results. 

Index Terms 

1^ Connectivity, percolation, outage, MIMO, diversity, power scaling. 
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I. Introduction 

Multihop relay networks have received a lot of attention recently due to their ability to improve coverage 
^ and, thus, capacity in a geographical sense. Many of these networks - such as mesh networks, vehicular 
^ networks, wireless sensor networks, and ad hoc networks - possess commonality insomuch as the number 
^ and distribution of nodes in the network is often random. A considerable amount of research on random 



O 



networks has been conducted in the past (see, e.g., [[T|-[|3|). From a communications perspective, it is 
of paramount importance to understand the connectivity properties of such networks. This understanding 
2 can lead to improved protocols and network deployment methodologies in practice [4|. 
CN In recent years, researchers have adopted and adapted a number of techniques from the physics and 
^ mathematical communities to study connectivity in random networks. Perhaps the most directly applicable 
. ^ theory that has been employed is that of continuum percolation [5], owing to its long and rich history of use 
^ in describing particle clustering in statistical physics and fluid dynamics [|6|-p0|. In particular, percolation 
theory is concerned with the emergence of a single large connected component (possibly in addition to 
other finite connected components) in a large (typically unbounded) graph, and dictates the minimum 
node density - which is known as the critical node density - required to obtain such a component. The 
links to network connectivity are obvious. Percolation theory has been applied in recent years to identify 
and analyze power management techniques that can be used to ensure network connectivity [11 1, p2| . 



S3 



It has also been used to demonstrate the benefits that node cooperation gives to improving connectivity 



[13|-[16|. Further applications of the general theory of percolation (including bond and site percolation) 



can be found in studies of network resilience p7| , hybrid networks (i.e., random networks with a regular 



element) [18], coverage and connectivity in wireless sensor networks [19[ (a direct application of [9[), 
and the information theoretic capacity of networks [ |20| . 

While the benefits of using percolation theory to explore asymptotic connectivity issues in random 
networks is clear, the theory does not directly address the question: what is the probability that all 
nodes in the network are connected? This question falls under the heading of full connectivity rather than 
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percolation. The answer to this question is, of course, related to a number of parameters, such as the fading 
environment, the path loss model, the node density, and the geometry in which the network resides. Many 
researchers have studied full connectivity, typically in some asymptotic regime. Two particular network 
models have been popularized for the study of full connectivity: the extended network model and the 
dense network model. The former relates to the case where the node density is finite and the network 
area is large, whereas the latter specifies a finite network area with a high node density (see, e.g., p7| 
and the references therein). 

With regard to extended networks, in [13|, p^, the authors assume a unit disk connection mo JeQ and 



give conditions on the path loss exponent for which full connectivity can be achieved in one-dimensional 
and two-dimensional networks. For dense networks, a number of scaling laws have been published. For 



example, in [122 1, the authors derive a power scaling law that ensures full connectivity is achieved almost 
surely as the number of nodes in the network tends to infinity. In [23], scaling laws are given for the 
number of nearest neighbors (i.e., connections per node) that are required to achieve full connectivity 
asymptotically in the number of nodes. Related results are given in [ [24 | for sectorized networks. 

More practical connection models (as opposed to the unit disk model) have also been considered in 



the literature [25[-[|30]. Specifically, the authors of [ |26| , [27 1 considered a probabilistic connection model 



whereby log-normal shadowing was incorporated into the path loss model. Both 1 -connectivity and k- 
connectivit)]^ were considered in those papers. In [29|, the authors considered the effects of fast fading (in 



particular, Rayleigh and Rician fading) on full connectivity in dense networks, and derived expressions 
for the node isolation probability for several cases of interest. Note that this probability effectively defines 
a first order approximation of the full-connectivity probability for dense networks, a point that we will 



elaborate upon later. Finally, the authors of [30[ have recently studied connectivity with respect to the 
decay properties of a general pair-connectedness probability function. 

The various contributions related to dense networks typically specify some geometry in which the 
network resides. These geometries are commonly taken to be squares or circles in two dimensions. One 
aspect that is common to virtually all of this research is that boundary effects are neglected in the interest 
of deriving a simple elegant result. In this paper, we show that such ejfects should not, in general, be 
neglected since they tend to dictate performance (in terms of the full-connectivity probability) in the high 
density limi^ We provide a constructive explanation of this assertion through the application of a novel 
cluster expansion model, first developed in [ ,3T| , [32|, which admits an accurate first order inhomogeneous 



approximation in the hmit of large density. Such models arise frequently in statistical physics to study 
the interaction between particles. The application of this approach to the problem of network connectivity 
carries many advantages, which together lead to the novel contributions of our work. These contributions 
can be summarized as follows: 

1) We utilize the inhomogeneous cluster expansion model to analyze the effects that boundaries have on 
connectivity for four important point-to-point link models: single-input single-output (SISO), single- 
input multiple-output (SIMO), multiple-input single-output (MISO), and multiple-input multiple- 
output (MIMO). To facilitate this analysis, we define the notion of the mass of connectivity, which 
is an average measure of pair connectedness in a given volume taking boundaries into account. 

2) We derive diversity and power scaling laws that dictate how boundary effects can be mitigated (to 
leading order) in confined dense networks. These laws are given for all four of the aforementioned 
point-to-point link models. 

3) We corroborate the general theory detailed in [ |3T| by analyzing boundary effects for dense net- 

'This model specifies that connection between two nodes is achieved if and only if the distance between them is at most r, which is some 
fixed positive number. 

fc-connected network is one that remains fully connected when k' < k point-to-point connections are broken. 
^See Fig. |l] for a qualitative illustration of this concept. 



Fig. 1. An illustration of the effects that the boundaries of the network domain have on connectivity. In this example, the network domain 
is a square of side length i = 10 in arbitrary units. The node densities are 0.5, 1, and 1.5, reading left to right. The positions of the nodes 
for each depicted realization are signified by pins, while the total probability of a new node introduced to the network connecting to any 
other node is indicated by the contours and shading (blue valleys indicate a low probability of connection while red peaks indicate a high 
probability). Notice that at high densities, low connection probabilities are concentrated near the comers. 



works comprising MIMO point-to-point links confined within a right prisrrQ This example is both 
instructive and insightful, illustrating the versatility and advantages gained by employing the new 
cluster expansion theory in network analysis problems. 
The rest of the paper is organized as follows. We begin in Section |II] by giving some preliminary details, 
which include the point-to-point path loss model and pair-connectedness probabilities for SISO, SIMO, 
MISO, and MIMO links. In Section III, we give details of the cluster expansion model and derive diversity 
and power scaling laws for mitigating boundary effects to leading order. We then go one step further by 
analyzing the connectivity of a network confined in a right prism. Analytical and numerical results for a 
specific right prism are given in Section IV-E Finally, we draw some conclusions in Section fV} 



II. Preliminaries 

In this section, we provide preliminary details of the pair connectedness model that we adopt. Specif- 
ically, we have that nodes i and j are directly connected with probability H(d(vi,Vj)), which we write 
as H (rjj) or just Hij, where the distance function d{ri, vj) is non-negative, zero only when rj = r^, and 
symmetric. We define the pair connectedness probability Hij as the complement of the information outage 
probability between nodes i and j. This well-understood metric is a natural choice for pair connectedness 
since it provides fundamental insight into network connectivity behavior, and can even be transformed 
directly into a unit disk model through appropriate parameter definition if so desired; however, it should 
also be noted that other pair connectedness models can easily be chosen, such as a model based on the 
average bit-error rate of a point-to-point link. The outage probability is generally parameterized by the 
received signal-to-noise ratio (SNR), which in turn is a function of many system parameters as well as the 
path loss. Consequently, we build our model of pair connectedness starting from a fundamental definition 
of the path loss model we employ. 



A. Path Loss Model and Outage Probability 

The received power of an electromagnetic wave decreases with distance like where r] is an 
environment-dependent decay parameter. Typically, r] = 2 if propagation occurs in free space, with r] > 2 
in cellular/cluttered environments or through walls (see, e.g., [33| and references therein). It follows that 
the SNR at the receiver (assuming a fixed transmit power Pt and a sufficiently narrow bandwidth]^ also 

''a right prism is a polyhedron constructed by taking an n-sided polygon as its base, replicating it and translating it in the "vertical" 
direction, then connecting the corresponding sides. Thus, right prisms are representative of, for example, many room configurations in a 
building. 

^Wideband channels sometimes exhibit different path loss behavior, which varies as a function of frequency. Discussion of these channels 
is beyond the scope of this paper; the interested reader is referred to I34j, (35) and references therein. 
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decays like r ^. Now, the outage probability for the link between nodes i and j is defined as the probability 
that the ij link cannot support a given rate Ro in bits per complex dimension, which can be written as 

Pi, = Ft (log2 (1 + SNR ■ Xij) < R,) 

where Xij denotes the random variable signifying the power of the channel between nodes i and j. The 
pair connectedness probability is simply Hij = 1 — Pjj, which, after some manipulation, yields 

(r) = 1 - Fx,, {f5r^) (1) 

where F^,^ is the cummulative distribution function of X^j, and /3 is a constant - which depends on the 
frequency of the transmission, the power of the noise process at the receiver, and the transmit power - 
that defines the length scale. It is important to note that (3 is inversely proportional to the average received 
SNR. Finally, we point out that by letting — )■ oo, we obtain the unit disk connection model: 



H,, {r) = { 1 - Fx,, , r = 1 . (2) 




B. Point-to-Point Link Models 

We consider four general point-to-point link models: SISO, SIMO, MISO, and MIMO links. For 
simplicity, we consider the case where individual channel fading distributions follow a Rayleigh model, 
and all channels are statistically independent. It follows that Xij has a standard exponential distribution 
in the SISO case and a chi-squared distribution with 2m degrees of freedom in the SIMO/MISO cases 
(where m is the number of diversity branches employed) p6J . Thus, for SISO links, we have 

H,, (r) = e-P^' (3) 

and for SIMO/MISO links, we have 

r (m, /3r^) 
r (m) 

where F (a, x) is the upper incomplete gamma function. Equation ^ reduces to the SISO case when 
m = 1. 

For MIMO links with rit transmit antennas and rir receive antennas, we are mostly concerned with 
ensuring connectivity is achieved. Thus, we assume beamforming is applied at the transmitter of each 



H^, (r) = ^^^7^ (4) 



node while maximum ratio combining (MRC) is employed at the receiver [37 1. For this so-called MIMO 
MRC channel, the pair connectedness probability becomes 

Hij (r) = 1 - Km,n det (7 (n - m + i + J - 1, (3r'^))^^ (5) 

where n = max(nt,nr), m = mm.{nt,nr), i,j = l,...,m, 7(a,x) is the lower incomplete gamma 
function, and 



-1 

{n — i + 1)T (m — i + 1) 



In order to aid analysis, we will focus on the special case where m = 2. This restriction is justified on 
the basis of pragmatism: indeed, one would envisage that low complexity may be a requirement of nodes 
operating in dense networks. This is certainly the case in wireless sensor networks where sensors are 
often powered by batteries ll3|. In any case, we maintain some level of generality by not restricting n. 
For m = 2 and general n, we can express Hij as 

Hij (r) = 1 - nP {n - 1, /3r'') P (n + 1, jSr'') + (n - 1) P (n, jSr^'f (6) 

where P (n, x) is the regularized lower incomplete gamma function. We will use these formulae for Hij 
in the next section to analyze the probability that a network is connected at the boundary of the confining 
geometry. 
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III. Probability of Full Connectivity 

In this section, we develop a novel theory of the probability of full network connectivity in confined 
geometries through the use of a cluster expansion technique. Cluster expansions are frequently used in 
statistical mechanics and fluid dynamics to study the interaction between particles. In our application 
of network connectivity, we first give an overview of the model, then provide details of the first order 
expansion, which can be used to study connectivity in the high density limit for general pair-connectedness 
models. For more details on this model, the interested reader is referred to p2| . We then utilize the pair- 
connectedness functions given in the previous section to analyze the probability of full connectivity at 
the boundary of the confining geometry to leading order, providing rules that dictate how diversity and 
power can be scaled to mitigate boundary effects. 



A. First Order Cluster Expansion 

Consider randomly distributed nodes with locations G V C M*^ for i = 1, 2, . . . , according to a 
uniform density p = N/V, where V = \V\ and |-| denotes the size of the set. Here, we use the Lebesgue 
measure of the appropriate dimension d. We define the average of a quantity as 

(A) = I A (ri, r2, . . . , r^) dridrs ■ ■ ■ drAr. 

We define some useful notation. Let S = {1, 2, 3, ... , A^}. A graph g = (A, L) consists of a set A C S* 
of nodes, together with a collection L C {(i, j) E A: i < j} of direct links, that is unordered distinct 
pairs of nodes. As a slight abuse of notation, we write (i, j) G (7 to denote that (i, j) is an element of the 
set of links L associated with the graph g. We write for the set of graphs with nodes in A, and G"^ 
for the set with nodes in A and largest connected component (cluster) of size j with 1 < j < \A\. 

The probability of two nodes being connected or not leads to the trivial identity 

1=H,, + {1-Hij). 

Multiplying over all links with nodes in a set A expresses the probability of all possible combinations. 
This can be written as 

1 = n + (1 - ^^^■)) = E ^5 

i,jeA;i<j geG^ 

where 

n n (^-H^J)■ (8) 

The sum in ^ contains 2l'^l^l'^l~^^/^ separate terms. Setting A = S, this can be expressed as collections 
of terms determined by their largest cluster, which yields 





P/c(ri,...,riv) a<,(l-^^^i) 



For a given configuration of node positions rj G V, assuming that the nodes are pairwise connected 
with independent probabilities Hij, the first term in equation (|9]) is the probability of being fully connected 
P/c(ri, . . . vn). The average of this quantity over all possible configurations Pjc = (P/c(ri, . . . r^)) is the 
overall probability of obtaining a fully connected network and is our desired quantity of interest. Hence, 
rearranging equation (|9]) allows us to obtain expressions for Pj^ in a consistent way while keeping track 
of correction terms. 

Note that in the high density limit of p — 00, the right hand side of (|9]) is dominated by the first term, 
which yields Pfc ~ 1, and hence the network is fully connected with probability one as expected. The 
approximation symbol is used here to indicate that first and higher order corrections are being ignored. 
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The first order approximation is obtained by expanding the second term on the right-hand side of Q 
explicitly, taking into account all possible ways of getting an — 1 cluster. Thus the average probability 
that an A^-node network confined in V is fully connected (to first order) is 



P 



1 - 



1 - 




N 





N 



[ n (l-/f(r,,v))dri---dr;v 



(10) 



where the simplification in the third equality is valid since all nodes are identical, meaning the sum over 
p can be factored out. 

At this point, we deviate from common practice and assume that the network is not translationally 
symmetric, i.e., the system is inhomongeneous. This ensures that we will include boundary effects in our 
analysis. Indeed, it turns out that it is not only important to include such effects, but that they dictate 
performance in the limit of high density. This is demonstrated to a large degree through the first order 



cluster expansion. In particular, we can progress from (10) to obtain 



P 



fc 



N 
V 



1 

V 



i^(ri2)dri 



N-l 



1-p f e-P^vH{^i2)dr^ (1 + (A^^^)) dr2 
Jv 



(11) 



in the limit of large N . In fact, this equation was recently given in pO[ (eq. (8)), where V was scaled 
exponentially with p, which effectively implies boundary effects are ignored. For our approximation, 
however, we only require that V ^ p, ov equivalently V ^ (cf. [32 1 for more details on scaling 
limits of our theory). This is a key difference from previous approaches reported in the literature. 

Equation ( [TT] ) suggests that in the high density limit, the probability of having a single — 1 connected 
cluster {i.e., having an isolated node) is dominated by nodes that are situated in "hard to connect" regions 



of the available domain V. This is because the outer integral in ( 11 1 is dominated where the integral in the 
exponential is small, which occurs at corners, edges, and faces. For nodes located near these geometric 
effects, the volume in range of the nodes is small. 



B. Pair Connectedness and Scaling Laws 



Taking a closer look at the integral in the exponent in ( 1 1 ), we see that it effectively defines the mass 



of connectivity for a given pair-connectedness function at Y2, and we label it accordingly: 

Mh{v2) = I H{ru)dr,. (12) 



V 



Indeed, Mh can be studied in detail for different pair-connectedness functions, and power and diversity 
scaling laws for bounded network connectivity can be gleaned from such an analysis. In the discussion 



7 



that follows, it will be understood that Mh is a function of V2, and thus we will refrain from explicitly 
writing out this relationship for the sake of brevity. 

We now incorporate the pair-connectedness functions given in Section [II] into the model. A common 
feature of these connection functions is that they decay exponentially with r. This allows us to separate the 



dri integral in (11) and extract its leading order behavior by supposing that r2 is located somewhere on 



the boundary dV of the network domain V C M [32|. Thus, for a general pair connectedness probability, 
we have to leading ordeij^ 

Mh=( r r^-^Hi, (r) dr\ ( [ do] (13) 



where = 27r'^/^/r (c//2) is the full solid angle in d dimensions, u = J d^l E (0,i7) is the solid 
angle as seen from r2, and we call Af^ the homogeneous mass of connectivity since it characterizes the 
mass of connectivity in a homogeneous system. Note that the domain of integration in the dr integral is 
not truncated from above because the integrand is exponentially decaying and, thus, does not contribute 
significantly to the integral for large arguments (i.e., for nodes located far from the node at r2). This 
approximation clearly demonstrates the influence that boundaries have on connectivity. Specifically, by 
writing 

Mh = SQM'h (14) 

where 5 = u/Vt E (0, 1), we see that connectivity is dominated from regions that are hard to connect and 
is controlled by the solid angle available to them. However, we cannot ignore the contribution of M^, 
which will be seen to be useful in determining ways to mitigate a small solid angle contribution. 

We now investigate the leading order behavior for the four link models detailed in Section |ll| placing 
particular emphasis on how diversity and power control can be used to mitigate boundary effects. Note 
that in order for us to obtain a more detailed view of boundary effects, we must expand ([T3|) beyond 



leading order; a discussion to this effect is given in Section IV 



1) SISO, SIMO, and MISO Link Models: For SISO, SIMO, and MISO pair connectedness, we have 

M^-fr-.^^dr (15) 

Jo r (m) 

The dr integral can be calculated directly by using the standard integral definition of the incomplete 
gamma function 

T(m,x)= / r-^e-Mt 



and exchanging the order of integration, which eventually yields 

M'h = ) (16) 

P^dV (m) 

For uncluttered (r/ = 2) d-dimensional networks with point-to-point SISO links, m = 1 and we have the 
simple result 

Intuitively, one might expect that boundary effects can be mitigated somewhat through the use of 



diversity (dictated by m in this case). Expanding (16) to leading order as m grows large yields 



d 



M'^ = ^{l + 0{m-')). (17) 

/3^d 



'We emphasize that equaUty is only to leading order here. 
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Fig. 2. An illustration of the M'^ scaling law for SIMO/MISO systems. The exact expression given by 1 16 1 is represented by solid lines 
for d = 3 and 77 = 2, 3, 4 moving from top to bottom. The dashed curves correspond to the leading order behavior. 



Thus, we see that although an increase in m results in better connectivity, the rate of increase is highly 
dependent upon the dimension of the space in which the network resides and the path loss exponent. 
Indeed, by scaling m, it is possible to obtain a progressive improvement in performance in uncluttered 
environments {rj ^ 2), whereas cluttered environments {rj > 3) will typically yield diminishing returns. To 
illustrate this scaling law, we have plotted the exact expression for M'^j along with the leading order term 
in Fig. [2} Not only is the leading order behavior apparent in this example, but we see that the leading 
order, in fact, provides an excellent approximation. 

As a final note, recall that (3 is inversely proportional to the average received SNR at a given node. 
Consequently, it can be seen from ( [T6] ) that an increase in the transmit power for each node has the same 
scaling effect as an increase in the diversity order m. However, it should be noted that this analysis is only 
valid when such an increase is less than exponential relative to the system size, i.e., scaling the transmit 
power must not counteract the exponential decay in pair connectedness. 

The fact that boundary effects can be somewhat mitigated through the appropriate use of diversity 
and/or power scaling is a powerful conclusion with significant practical implications, particularly since 
these scaling laws have identical order. For example, some networks may be energy/power constrained 
(e.g., wireless sensor networks and cognitive networks) whereas others may be constrained by form size 
(e.g., some vehicular networks). The scaling laws presented here provide engineering insight into how to 
design networks for a plethora of scenarios. 

2) MIMO Link Models: We can perform the same analysis for MIMO systems. In this case, we restrict 
our attention to systems where m = 2, for which Hij can be rewritten from ([6]) into the form 



H^, (r) 



nV (n - 1, /3r^) - 2T (n, (3r'^) + (n - 1)"^ T {n + 1, (3r'^) 



rfn- 11 



r {n, /3rY - r (n - 1, /3r^) T{n + 1, /3r^) 

^ r (n) r (n - 1) • ^ ^ 
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Using this form of Hij, we can evaluate by integrating term by term and using recurrence relations 
of hypergeometric functions to yield 



M',= - "i ^ ^ + -^4 y-r-F(n-l^2n+'--n + l--l 



V 



n — 1 ^ ( d d d , , ^ 

F n-l + -,2n+-;n + l + -;-l (19) 



where F (a, b; c; z) denotes the hypergeometric function. The terms involving hypergeometric functions 
result from integrating the quadratic terms in the expression for Hij given above. Again, the effect 
of the solid angle defining the volume in range is apparent. Although the expression given above is 
quite cumbersome and does not necessarily lead to much intuitive insight, we can easily evaluate the 
approximation for specific values of n. For example, we have for n = 2 

= ^ , ^—^ . (20) 

On the other hand, we can obtain a scaling law in n by noting that for large n, the contributions of the 
P-functions in (j6]) are not significant for (3r''^ < n. Moreover, for (3r''^ > n, Hij (r) decays exponentially. 
Thus, we can approximate Hij (r) by a step function with a transition at r = (n//?)''. We present an 
argument for making this approximation in Appendix lAl It follows that we can write 



1 







d \/3 



d 



M'h= r' r''-'dr + e{n) = ]('^Y +e{n) (21) 



where the error term is given by 



n \ rj 



1 

OO 







e (n) = / r'^-i {Hij (r) - 1) dr + / , r'^-^Hij (r) dr . (22) 



^ 



Note that e_ (n) corresponds to the negative contribution and e_|_ (n) is the positive contribution. It remains 
to determine the order (in n) of the error term. In fact, it can be shown (by way of accurate approximation 
if not rigorously) that 

e{n) =0 (nn-'A . (23) 



A derivation of this result is given in Appendix |B] Thus, we can write 



M'^ = ^(l + 0(n-'^)) (24) 

I3^d ^ ^ 



which is completely analogous to the SIMO/MISO case. This scaling law is illustrated in Fig. |3} where 



we have plotted the exact expression for given by (19) along with the leading order term. The leading 
order growth in n is apparent in this example, but in contrast to the SIMO/MISO case, it is not a good 
approximation. This results from the fact that the first correction term is O (n^^/^) rather than O (n^^) as 
is the case in SIMO/MISO systems. Finally, it is easy to see (by observing ( [T9] )) that the same transmit 
power scaling law holds, i.e., M'^ scales like (3^^ . 
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Fig. 3. An illustration of the M'u scaling law for MIMO systems. The exact expression given by \\9\ is represented by solid lines for 
d = 3 and t) = 2, 3, 4 moving from top to bottom. The dashed curves correspond to the leading order behavior. 



IV. General Formula and Applications to Right Prisms 

In fy2)\, we developed a general formula for studying the full-connectivity probability in the limit of 
high density for any network domain. This formula takes the form 

d 

i=Q ji 

where, again, d is the dimension, Gj^ is a geometrical factor for each object ji of codimension i, and 
Vj. is the corresponding d — i dimensional volume of the object with solid angle coj^. In the spirit of 



this formula, we now derive an expression for Pf^ by evaluating (11) when the network in question is 
located in a right prism and the point-to-point links are modeled as uncluttered (i.e., i] = 2) 2 x 2 MEMO 
channels. The choice of MIMO point-to-point links for this analysis was taken to demonstrate that our 
methodology can be applied to reasonably complicated scenarios. Indeed, SISO, SIMO and MISO links 
can also be analyzed. 

The pair connectedness probability for 2x2 MIMO channels with 7] = 2 simplifies to 

H (r) = e-^'" ( /3V" + 2 - e~^''') (26) 



where r is the distance between two nodes. To derive P/c in the form of (25), we evaluate the integral 
(cf. ^) 

I e~''^V^^(ri2)dri(Jj.2 ^ I g-Mfe(r2)jj.2 (27) 

Jv Jv 

for each local feature of the geometrjj^ This illustrative example will further show how geometric effects 
influence connectivity. Moreover, the choice of a right prism as the confining geometry demonstrates the 



^Here, we have explicitly indicated that the mass of connectivity Mh is dependent upon r2 since we will be concerned with first order 
correction terms. 
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power and versatility of our theory, particularly since many practical geometries can be well approximated 
by such a polyhedron. 

At this point, it is beneficial and instructive to give an outline of the approach we use to analyze 
networks confined in right prisms. Our general method is to begin by considering features with the lowest 
dimension, i.e., comers in this case. We then move to edges, then faces, and finally the bulk of the prism. 
At each step, we ignore objects of lower dimension since we have accounted for their contribution to 
connectivity in previous steps. It will be observed that this is a particularly powerful approach when 
analyzing the effects that the bulk and faces have since the surface (volume) of a right prism is locally 
equivalent to that of a sphere of the appropriate surface area (volume). We now give details of our analysis, 
and follow this discussion with a specific example. 



A. Corners 

For a right prism, each corner is the product of three intersecting planes that can be oriented such that 
at least one edge connected to its vertex is normal to an adjoining face. This suggests that we should use 
cylindrical coordinates to perform calculations for this case, where we locate the origin at the vertex and 
the 2;-axis is oriented along the edge that connects the two identical polygons. 

The distance between two points in cylindrical coordinates is given by 



d (ri, r2) = yrl + rl- 2rir2 cos (^^i - 62) + (21 - 
where (rj, 6'j, zi) are the coordinates of a node located at rj. Thus, in order to evaluate the inner integral 



of (27) near the comer, we let r2 be located near the corner and expand H {^12) near r2 = and Z2 = 0, 



which, to first order, yields 



1 



H (ri2) = e-'^H+^?) (e^('-?+^?) (2 + (3^ (r? + zl)' 

+ 2/3e-2^('-'+^') [2 + (3 [rl + zl) {(3 [rl + zl) - 2)) - 2^ 

X {Z1Z2 + rir2 cos {61-62)) . (28) 



We can now perform the inner integrals in ( [27] ) as follows: 
MHiv2)= / / / nHiv,2)dnd6^dz, 




+ Tvrrs (sin ^2 - sin (^2 - ^) ) ) (29) 

where tD is the angle of the corner with < < vr. Note that semi-infinite integration is allowed here 
due to the fact that H is exponentially decreasing but the system size is largej^ Thus, contributions to the 



integral at distant boundaries are negligible. Using (29), we are now in a position to calculate the outer 



integrals of (|27j), which yields 

g-pM^(r2)^j.^ = / / / r2e-^^V^(ri2)dri^^2^^2^^2 




256/3^ CSC ^9 _(^^z^!V^P^ 
= e 16/3^/2 (30) 

The regions of integration are the same here as for the inner integrals; however, note that the order of 
integration changes. 

All that remains is to enumerate the 2n comers for a prism constructed from an n-sided polygon. Each 
comer is defined by the angle d. For example, a cuboid, which is a right prism formed by replicating and 
translating a square, has eight comers, all of which have angle d = tt/2. 

*In particular, if L denotes the typical length of the geometry, then we require \/T3L 3> 1. 
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B. Edges 

Now we consider geometric features of dimension one: edges. Let L be the length of the edge in 
question. The calculations for this case are also facilitated by the use of cylindrical coordinates, but where 
the origin is located at the center of the edge. Thus, the corners are located at ±L/2 and the angle of the 
comer is Since we wish to ignore effects from comers, faces, and the bulk, we expand H about r2 = 
and Z2 = 0, which gives (28). Calculating the inner integrals in (27) to first order yields an expression 

that has exp (— /3L^/4) and erf {^L^J (3/2^ terms. Again, assuming that ^/]5L ^ 1, we can make the 

approximations exp {—(3L'^ /A) ^ and erf (^L^/f3/2^ ^ 1. It follows that Mh (r2) can be evaluated to 
yield 

Mh {r2) = j- I ^^^-^./^^ + 77rr2(sin^^2-sin(02-^)) 1 (31) 




/3 

where the integrals are performed over ri G (0, oo), 6i G (0, ■&), and zi G (— L/2, L/2). The outer integrals 



in (27) can now be performed to yield 



Again, all that remains is to enumerate the 3n edges. 



C. Faces 

For the contribution of the faces to the full-connectivity probability, we employ a local equivalence 
argument that allows us to greatly simplify the analysis. Specifically, we have covered the corner and edge 
calculations above, and thus we ignore contributions from these features when considering faces. Thus, 
one can imagine deforming a prism of surface area S into a sphere of the same surface area, the radius 
R of which is defined by the relation S = AnR^. For a general right prism, the surface area is given by 

S = 2B + ph 

where B is the area of the base {i.e., the n-sided polygon), p is the base perimeter, and h is the height. 
If the base is a regular n- sided polygon with side length s, we have 

b = —s cot — h nsh. 
2 n 

Thus, this argument allows us to treat any convex right prism that we wishj^ 

Using spherical coordinates along with the fact that the distance between nodes at ri and Y2 is given 

by 



d (ri, r2) = \l ^1+ 2rir2 cos 6 

where = |rj| and 6 G [0, vr] is the angle between the nodes, we expand H near the surface of the sphere 
{i.e., r2 = R) and perform the inner integrals in (27) to obtain 

AfH(r2) = 2jr/ / if sm9J/(r,2)dMri 




'%14(fl-r,)| (33) 



'Convexity is required since we only consider line-of-sight connections between nodes, although small-scale scattering is accounted for 
through the chosen pair-connectedness models. 
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where the factor of 2n follows from integration over the azimuthal angle. To arrive at this result, it was 
assumed that a//?-R ^ 1, which allows us to make similar approximations to the error functions of the 
form erf {cy/]3R) and exponentials of the form cxp {—c(5R^) for some constant c > as was done for 
edges. Furthermore, we have made the additional approximations c/R and c/R^ ^ 0. The outer 
integrals can now be performed to yield (to dominant term in R and p) 

j-R rn ^RR'^ (23-V2)7r3/2p 

27r / / r^sin^ e-^^^^'-^MMrs^ ^V^e . (34) 

Jo Jo iP 

Generalizing this result to any right prism, i.e., substituting S = AnR'^, gives 

/ e-^^«("-2)dr2 = — e ^^"^ . (35) 
Jv 77rp 

The faces do not need to be enumerated in this case since we have accounted for all faces through the 
local equivalence argument. 

Finally, it is worth mentioning that face contributions can also be calculated in a lengthy manner by 
using cartesian coordinates. This works well for rectangular sides; however, one must be more careful 
when considering more general n-sided polygons. In any case, it is straightforward to show for certain 
simple cases that the proposed approach to the calculation yields identical results to the more complex 
cartesian approach. 



D. Bulk 

For the bulk contribution, we apply the same local equivalence argument that was used for the face 
contributions, but where the expansion in H is performed at r2 = 0. In other words, we consider a sphere 
of radius R determined by the relation V = |7ri?^, where 

V^Bh 

is the volume of the right prism containing the network. For a prism constructed from a base that is a 
regular n-sided polygon with side length s, we have 

V ^ -hs^ cot-. 
4 n 

Expanding H about the origin to first order, we can perform the inner integration to obtain 



poo pTT 

Jo Jo 
_ {23 -V2) /tt^ " 



4 

Finally, the outer integrals can be evaluated to yield 

-R PTT (23-^/2)^3/2p 



(36) 



27r / / rlsine e-P^"^'''Medr2^Ve . (37) 

Jo Jo 
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E. Example: the "House" Prism 

We now apply the calculations detailed above to a particular scenario. Specifically, we consider a 
network comprised of 2 x 2 MIMO links with rj = 2 operating in a right prism that resembles a typical 
"house", as illustrated in Fig. |4j The base of the house is a square of side L. The total height of the 
prism is 3L/2, and the apex is right angled. This example somewhat illustrates the versatility of our 
theory, particularly since most previously published results have considered very simple geometries, such 
as squares and circles in two dimensions, or cubes and spheres in three dimensions. It should be noted 
that this is an arbitrary example, and in fact many more complicated geometries can be analyzed with the 
proposed methodology. 

Analysis of the "house" prism is straightforward using the calculations given in the preceding sections. 
We first identify the different comer angles and count the multiplicities of each. There are ten corners in 
this prism, six of which have angle ■& = j2 and four of which have angle = Sir /A. Thus, the following 
two contributions to the general formula for Pj^ arise from the comers (cf. (30)): 

512/33 (23-^.^x3/2 

and 

1024\/2/3^ {23-V2)3p (^^3/2 

Next, we perform a similar step for the edges, of which there are fifteen. Thirteen edges are right 
angled: nine of these have length L while the remaining four have length L/^/2. The other two edges 
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TABLE I 

Contributions of the various geometrical features of the "house" prism to the general formula for Pfc given by 



(25 (. The angle j9 = n/2 for type-1 corner/edge contributions, and i9 = 37r/4 FOR TYPE-2 contributions. 



General Formula Parameter 


Corners 


Edges 


Faces 


Bulk 


Volume (VjJ 
Solid Angle (tOj^) 
Geometrical Factor (Gj- ) 


1 

■a 

256/3^ csciJ 


r ^ 
2& 

16/3^ CSC i5 
497r^ 


S 
2tt 

2/3 
Ttt 


V 
4ir 
1 



have angle ?9 = 37r/4 and length L. Thus, we can write the following two edge contributions to the high 



density expression for Pfc (cf. (32)): 



and 



E, = L{9 + 2V2)^^e-^^i^) (40) 

16^2/32 (23-^)3p .^x3/2 

E2 = 2L , , , e 55 — UJ . (41) 



497r2p2 

For the faces of the prism, we calculate the total surface area to be 

ll + 2v^^2 



S = ^r^L\ (42) 



Thus, using the local equivalence argument presented above, we can substitute this area into (35) to obtain 
the contribution of all of the faces to Pfc. We denote this contribution as F. Similarly, we can take the 
volume of the prism to be 

V = ^L^. (43) 

4 

Substituting into (37) yields the bulk contribution to Pfc, which we denote by U. 



Finally, the general formula is obtained through the summation of all contributions, which yields 

Pfc^l- p{Ci + C2 + Ei + E2 + F + U). (44) 



It can easily be seen that this formula has the same structure as (25). In fact, the various contributions to 
the general formula are outlined in Table |I} 

Letting (3 = 1 and L = 7 for simplicity, we have plotted the general formula for this example along 
with numerical results obtained through simulations in Fig. [5} Furthermore, to illustrate the accuracy of 
our results at high density, we have plotted the probability of network outage (i.e.. Pout = 1 — Pfc) 
in Fig. [6j In these figures, the two solid curves (reading left to right in both figures) are the numerical 
and full analytical results (including contributions from all geometrical boundary features). Note that good 
agreement is achieved at high densities as expected. The left-most dashed curve (red) is the general formula 
with only a bulk term (i.e., all boundaries are ignored). This curve representes "conventional wisdom", 
where boundary effects are neglected, and serves as a benchmark for our results. This benchmark is overly 
optimistic and should be treated with caution. This is particularly clear in Fig. |6j The next dashed curve 
(yellow) includes bulk and face contributions. We see that by including even two-dimensional boundary 
effects, the model is improved significantly. The right-most dashed curve (green) includes bulk, face, and 
edge contributions, but not corners. 

Although these results show that our theory is in excellent agreement with simulations for high densities, 
they also illustrate the limitations of the theory for low to mid-range densities. However, one notices that, 
for example, boundaries are not as significant at low densities. Moreover, for mid-range densities, one 
may only need to consider boundary effects due to faces and/or edges, but not corners, in order to obtain 
accurate connectivity predictions. This observation can be explained in, perhaps, two ways. Firstly, from 
a mathematical point of view, our analysis relied heavily on expansions, both near boundaries and at high 
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Fig. 5. Analytical and numerical results for the connectivity of a network in a typical "liouse". The two solid curves (reading left to right) 
are the numerical and full analytical results. The left-most dashed curve (red) is the general formula with only a bulk term (i.e., face, edge, 
and corner terms are ignored). The next dashed curve (yellow) includes bulk and face contributions. The right-most dashed curve (green) 
includes bulk, face, and edge contributions, but not corners. 
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density. To maintain simplicity and tractability, these expansions were typically only to first order. This 
may be a limitation when considering low to mid-range densities. Secondly, one may intuitively expect 
that away from the region of high density there is a low probability of a node being located very close to a 
boundary, which of course is where our analysis is valid (again, expansions were taken near boundaries). 
In any case, our results actually suggest that the theory can be modified somewhat to be valid in regions 
of low and mid-range densities. This is certainly of interest, and is currently ongoing work. 



V. Conclusions 

In this paper, we analyzed the probability of full connectivity for dense networks confined within certain 
geometries. In contrast to much of the results reported in the literature, we showed that boundary effects 
should not be ignored, and in fact dictate performance in the high density regime. We built upon our cluster 



expansion model of full network connectivity, first reported in [31 1, by adopting four important point-to- 
point link models for pair connectedness: SISO, SIMO, MISO, and MIMO. This approach allowed us to 
study the effects that the boundaries of the confining geometry have on practical network connectivity, and 
led to the derivation of diversity and power scaling laws that dictate how these effects can be mitigated. 
Finally, we provided a practical example by analyzing the full-connectivity probability of a network 
comprising MIMO point-to-point links confined within a right prism. This example demonstrated the 
versatility of our theory in network analysis problems, and we hope that this theory will inspire many 
researchers in their own problems related to network design and performance analysis. 



Appendix A 

A.RGUMENT FOR APPRC 

Let H be given by 



Argument for Approximating Eij by a Step Function 



U ix) = 1 — nP in — 1, x) P (n + 1, x) + (ra — 1) P (n, x) 



First, consider the case where x ^ n. Employing the asymptotic equivalence [38, 8.11.5] 



P (n, X) ~ , . 

we can write 



Hix) r^i-- (: 



27rn Vn/ 

and thus (x) — t- 1 as — )■ cxo. 

Now, consider the case where x is close to but less than n. Let x = \n < n for some A ~ L Using 
the asymptotic relations [^38^ 5.11.3, 8.11.6], we can write 

7(?T,,x) {n — x)~^ x"'e~^ (Ae^'^^'^)" 



rH y2^r^"~5e-" v/2^(l-A)' 



P {n, x) 

But < Ae^-^ < 1 for < A < 1, and thus P {n, An) for fixed A < 1. It follows that H (x) ^ 1 for 
X < n and n large. 

Finally, consider the case where x is close to but greater than n. Let x = \n > n for some A ~ L 
Now we employ the asymptotic relations [j38| 5.11.3, 8.11.7], which give 



. , r(n,x) (Xe'-^y 



r(n) v^(A-l) 

But < Ae^"^ < 1 for A > 1, and thus P {n, Xn) 1 for fixed A > L It follows that if (x) for 
X > n as n — 7- oo. Thus, it is logical to make the stated step function approximation with a transition at 
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Appendix B 
Order of the Error Term 

To determine the order of e (n), we will make use of the expansions (3S^, 5.11.3] 



and [38 8.11.12] 



r (n) ~ y2^n"-^e-" (l + O (n"^)) 



(45) 



(46) 



Moreover, we have the following identity (which holds for any a, but where the logic behind the proof 
is given for a > 0): 

r {n + a,n) T [n + a, n)r [n + a) 



Tin) 



r{n + a) r (n) 
T (n + a,n + a) 



+ 



r (n + a) T {n + a) 



n+a 



r + a) 



Similarly, we have 



r 2n + a,2n 22("-i)+« i / ^ / 



(47) 
(48) 



We begin by determining the order of 



e+ [n) 



OO -| 



2l 1 T] 

PI 



if (a;) dx 



where H (x) = H ^(a;//3)^^''j . We can rewrite H (x) to yield the following bound on the pair connect- 
edness function: 

{x — n) x"'~^e~^ 



H (x) = 2Q (n, x)-Q {n, xf + - — + (1 - Q (n, x)) 



T(n) 



T{n) 



< 2Q {n,x) + 



Tin] 



T(n) 



-, X > n 



V 



where Q {n, x) = 1 — P (n, x) = T {n, x) /T (n). Now we can write 
1 



X'' Hh{x)dx 



— 

rjP-^V [n) 



d 



2i-2n---p i2n-l + -,2n 



d 



d 



— r {n) I (n, ra) + riF I n — 1 H — , n — ( 1 + 2- ) F I n H — , 



Using the expansions given in (45)-(48), it can be shown that this bound behaves like O [n^^ 2 
Considering the negative errorierm, we have 



e_ [n] 



— J / x"^Hc{x)<lx 
r]f3^ Jo 
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where Hc{x) = 1 — H (x), which behaves like o (x^" ^) for small 5 > as x — )■ 0, and increases sharply 
near x = n for large n (see Appendix [a]). Thus, we propose a semi-rigorous bound on He (x) in the form 
of an exponential with a slope at x = n matched to that of He (x), i.e., 

where a„ and 6„ are constants that may depend on n. Note that H^ {x) increases with x and, thus, has a 
maximum in the interval x G [0, n] at x = n. Furthermore, He {n) decreases with n, so He {x) is upper 
bounded by 

Furthermore, the slope of He (x) at x = n is given by 

^ng-2n (^2^" + e"7 + 1, n)) 



H'An) 



T{n + lf 



Matching the slope of the bound at x = n and setting the bound equal to He {n = 2) at the same point 
leads to the definitions a„ = He{n = 2) and 6„ = a^^H'e (n). 
Making the substitution t = n — x allows us to write 



— — 1 / \— — 1 



By the mean value theorem, we have 



n 



n)l<^^^^ / (1--1 e-'"'dt 



-1 d 1 

- — 1 /"n 



^1 d 



-1 



7 (1 - e-'"^) 

--1 

8^2^ (cosh 2 - 3)^ d_i f ^( -1 



for some G (0, 1). Thus, for d > r], |e_ (n)| = O (ra^ 2 
For d < r], we have 



d 1 



< -d + / x^e^dx. 

The first term on the right hand side of this inequality goes to zero exponentially with increasing n. 
Consequently, we focus on the second term. Notice that /„ (x) is convex, with a minimum at x = 

^1 — /hn- Moreover, /„ (1) — > 1 as n — )► 00. Thus, the major contribution to the integral /"/n occurs 

near the upper end of the integration region, i.e., x ~ To make progress, let y (x) = 6„x —{l — ^ log x 
such that /„ (x) = e^*^^^ and expand y to first order about x = n. This yields the accurate approximation 

In (x) ~ e3/(")+J''(")(^-") 
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We can now integrate easily to obtain 



-\ d d , , 

i — / 1 a 



&nri + J - 1 



and thus, to a good approximation, we can write 



\e. {n)\ < T + . ( - e'^-^"-)- 

2 



8V27r (cosh 2 - 3) 

~ ^ 1 '"''^ ^ 

and so we have shown that e (n) = O I n'' ^ 
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